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ABSTRACT 

A potential ens trophy conserving scheme and a modified 
a-coordinate have been coded in the UCLA general circu- 
lation model. The new model is being tested with global 
simulation experiments. 

Introduction 

During the past year we have modified several aspects of 
the UCLA general circulation model. Major changes are: 

The inclusion of a potential enstrophy conserving advec- 

tion scheme in the equation of motion; 

reformulation of the vertical differencing based on a 

modified a-coordinate in which the predicted PBL top is a co- 
ordinate surface; 

the inclusion of the turbulent entrainment parameteriza- 
tion reported by Randall at the Third NASA Weather and Climate 
Program Review. 

Other changes include: a potential temperature conserving 

horizontal advection scheme, a modified upstream scheme for 
advection of moisture, and inclusion of the ice phase in the cloud 
model used in the cumulus convection n?irameterlzatlon. 

The potential enstrophy conserving scheme 

Large truncation errors produced by poorly resolved steep 
topography can in a non-linear system affect large-scale waves. 

The convergence of the solution as the grid size is reduced will 
be slow unless the numerical scheme can prevent the false cascade 
of energy to small scales. In a purely two dimensional flow, 
governed by conservation of absolute vorticity, r)j ^ false 
cascade ma y be prevented by using a scheme that conserves 

enstrophy, , and kinetic energy. In the presence of topography 
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a two-dimensional bar o tropic flow is governed by conservation of 
absolute potential vorticity, q = n/7T, where tt is the surface 
pressure. In this case a numerical scheme that conserves 
potential enstrophy and total energy will guarantee an upper 
bound on the absolute enstrophy and once again prevent the 
cascade to small scales. 

As reported by Arakawa at last year’s NASA Weather and 
Climate Program Review, a preliminary study using shallow water 
equations with steep topography showed that a potential enstrophy 
conserving scheme was superior to a conventional scheme when the 
grid size was comparable to the scale of variation of the to- 
pography. These ideas have been extended to the general circula- 
tion model equations in the folowing way: The equation of motion 

and continuity in a 0-coordinate may be written in the form 

*1“ =“ r) k X V - 0 - V(<j> + h V»V) - 0aV7T 

g= - V.(Trf) -|^ (TO) 
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These equations guarantee conservation of i$q ttAo = ^ /ttA 0 for a 
layer between two material surfaces separated by an amount A0, 
provided V is independent of 0 within the layer, Va x ViT *= 0, and 
the vorticity is defined in terms of gradients at constant 0. 

The numerical scheme guarantees conservation under these 
constraints when the material surfaces coincide with its 0- 
coordinate surfaces. 

The modified 0-coordlnate system 

A major deficiency noted in the previous version of the 
UCLA general circulation model was a lack of organization in 
fields strongly influenced by the planetary boundary layer. 
Precipitation in the tropics, for example, was less sensitive to 
sea surface temperature than in either nature or earlier versions 
of the model. This was thought to be due to the indirect way in 
which the general circulation and the parameterized PBL were 
coupled. 


In the old model layers were defined by conventional 0- 
coordinate surfaces, as indicated by the heavy lines in the 
upper half of figure 1. The well-mixed planetary boundary 
layer, whose depth was predicted, was then allowed to occupy any 
fraction of the lower atmosphere. Typically, it was confined 
to a part of the lowest general circulation layer, as shown on 
the right of the figure, but it could also penetrate to higher 
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layers, as shown on the left. Coupling between this PEL and the 
general circulation model was handled by writing prognostic 
equations for the magnitude of the discontinuities at the top of 
the PBL and diagnostic relations for PBL properties. 

We felt that the difficulties arising from the weak 
coupling between PBL and general circulation model in this scheme 
may be overcome by making the top of the PBL a coordinate surface. 
We therefore redefined the vertical coordinate as follows: 


a 
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where p^ is a constant pressure above which we use pressure co- 
ordinates, Pg is the pressure at the PBL top, and p is the 
surface pressure. The PBL top then corresponds to a = 1 and the 
earth’s surface to a = 2. This coordinate allows us to make the 
PBL the lowest GCM layer, and PBL quantities become the prognostic 
variables from which the jumps at the PBL top are defined. To 
improve the accuracy of the Jumps, vertical resolution is in- 
creased immediately above the PBL. The model is being run with 
the 9-level (15-level when the stratosphere is included) structure 
shown in figure 2. 

These changes have been coded into the new version of the 
UCLA general circulation model. The model has been tested for an 
idealized limited domain and is now being tested with global 
simulation experiments. 
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OLD MODEL 



NEW MODEL 


Figure 1 - Schematic of the model's layers near the ground in the 
old and new 0-coordinates. Note that in the new version the PBL 
coincides with the lowest layer 
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